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ABSTRACT 



The solar tachocline, located at the interface between the 
latitude-dependent rotation of the convection zone and 
the rigid radiative interior, presents high gradients of an- 
gular velocity which are of particular interest for the mod- 
els of the solar dynamo and angular momentum trans- 
port. Furthermore, latitudinal and temporal variations of 
the tachocline parameters, if any, are also of particular 
interest in order to constrain models. We present a re- 
view of some of the theories of the tachocline and their 
predictions that may be tested by helioseismology. We 
describethe methods for inferring the tachocline parame- 
ters from observations and the associated difficulties. A 
review of results previously obtained is given and an anal- 
ysis of the new 6 years database of LOWL observations 
is presented which yields no compelling evidence of vari- 
ations or general trend of the tachocline parameters dur- 
ing the ascending phase of the current solar cycle (1994- 
2000). 
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1. INTRODUCTION 



The solar tachocline, so called after Spiegel & Zahn 
(1992), is defined as the layer at theTase of the convec- 
tion zone (hereafter, CZ) where important radial shear oc- 
curs. It is the transition zone where the angular velocity 
changes from its latitude-dependent value in the CZ to its 
constant and intermediate value in the radiative interior. 
We can emphasize three main reasons why this layer is of 
particular interest. 

(i)Shear turbulence and/or meridional circulation inside 
the tachocline may provide a mechanism for mixing ma- 
terial betw een the CZ and the radiative interior (e.g. B run 
et al., 1999; [Schatzman et al., 200CI| ) which is needed to 
understand, for instance, the burning of lithium and the 
depletion of helium and to reach a better agreement be- 
tween solar models and helioseismic observations (e.g. 
Richard et al., 19961; brun et al., 19991). 



rigid rotation rate of the radiative interior Hydrody- 
namical tra nsport by unstable shear flow (e.g. Chaboyer 
et al., 1995) an d transport by inter nal gravity waves in 
the t achocline ([Kumar et al., 1999|; K im & MacGregor, 
2000) have been studied but they have not been found 
efficient enough to lead to an uniform internal rotation 
at the present age of the Sun. This probably requires 
also the pres ence of an internal magnetic field (M estel & 
Weiss, 1987; |Charbonneau & MacGregor, 1993|; G ough 
&McIntyre, 1998). 

(iii) Finally, the tachocline is the best location for an os- 
cillatory solar dynamo which is generaly believed to be 
responsible for the solar magnetic cycle for the following 
reasons: (1) Its radial and latitudinal differential rotation 
has the ability to produce a toroidal field by shearing a 
pre-existing poloidal field. (2) The a-effect, the essen- 
tial mechanism for producing poloidal field from toroidal 
field, is usually located in the CZ but the tachocline can 
also prod uce a strong a-effect by magnetic buoyancy in- 
stability ( l^erriz-Mas et al., 1994 ) and/or by the unstable 
shallow-water modes ( pikpati & Oilman, 2000^ ). (3) Be- 
cause the tachocline (or part of it) may also be located in 
the slight sub-adiabatic overshoot layer, the toroidal fields 
can be stored for an extended period of time, and there- 
fore can be amplified and acted on by the a-effect before 
they escape to the surface through buoyant rise or be dis- 
rupted completely by convective shredding. 

We will first (Sect. ||) give some terminology about the 
different layers and their properties at the base of the CZ 
(also to be defined). Then, in Sect. ||, we present some 
models of the tachocline and their predictions that may 
be tested by helioseismology using methods presented 
in Sect. ^. Previous results obtained are summarized in 
Sect ||and we then use the new LOWL data to investigate 
latitudinal and temporal variations of the tachocline pa- 
rameters during the ascending phase of the current solar 
cycle (Sect. 0). 



2. LAYERS AT THE BASE OF THE CONVECTION 
ZONE: TERMINOLOGY, DEFINITIONS AND 
OBSERVATIONS 



(ii)The tachocline may be the seat for angular momen- 
tum transport processes that could lead to the observed 



As mentioned above, there is, beside the tachocline, an- 
other layer defined at the interface between the super- 
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criterion. In that case the downward convective energy 
transport {Fc < 0) occurs only in the other (lower) part 



F=F, 



Figure L Schematic representation of the overshoot 
layer between the radiative zone boundary Vrz ond the 
Schwarzschild boundary rsd where the temperature 
gradient equals both its adiabatic and radiative value 
V = Va = Vr j. Tcz represents the seismically observed 
base of the CZ. Two different models of the convective 
( Fc) and radiative ( Fr) flux variations are drawn for 
comparison. The total flux Ft — Fr. + Fr remains con- 



stant. The model (1) of Canuto (1997) allows a positive 
convective flux in the slightly sub-adiabatic and seismi- 
cally observed part of the overshoot in contradiction with 
the S chwarzschild criterion while the model (2) ofZ ahn 
(1991 ) is such that Fc < in all the overshoot layer 



adiabatic CZ and the radiative interior: the overshoot 
layer. Different models exist for this layer. The first 
particular radius, Vsd, corresponds to the Schwarzschild 
boundary which defines the beginning of the convectively 
unstable zone where the temperature gradient equal its 
adiabatic value V = Vad- A second particular radius 
is where the temperature gradient is taking its radiative 
value (V = Vr) and the total flux is equal to the radia- 
tive flux (Ft — Fr). This radius, r^z, defines the radiative 
zone boundary. The region between r^d and r^^ is what 
we here call the overshoot region (see Fig. |l|). The 'base 
of the convection zone' Vcz, measured by helioseismol- 
ogy, corresponds to the location where the temperature 
gradient changes abruptly. If we allow a slightly sub- 
adiabatic zone below Vgd then helioseismology cannot 
distinguish between the slightly super-adiabatic CZ and 
the slightly sub-adia batic overshoot layer (C hristensen- 
Dalsgaard et al., 1995; iGilman, 2000D and Trz <Tcz <1"sd- 
Calibrations made by pjasu (1997| ) lead to a value of 
Tcz — 0.713 ± O.OOli?0. Helioseismology can also place 
an upper limit for the slightly sub-adiabatic overshoot 
zone [rcz,rsd]- This upper limit is 0.05 pressure scale 
heig hts, about 2800Km or ^ 0.004^^ ( Basu & Antia, 
1997). This is therefore a very thin layer. Another ques- 
tion concerning this layer is about the sign of the con- 
vective flux inside. If we follow the Schwarzschild crite- 



rion, we should have i^c < as in Zahn (1991 )'s model 
because th is layer is (slightly) sub-adiabatic, but F erriz- 
Mas (1996) and [Canuto (1997D pointed out that, with a 
nonlocal mixing-length treatment of the CZ, this is not 
necessarily valid in the overshoot zone and we may well 
have Fc > and V < Vad in spite of the Schwarzschild 



of the overshoot zone I 



, called the ov ershoot layer 



proper ( I^anuto, 1997 ; Monteiro et al., 1998 ). In the same 
way, the zone with pos itive convective flu x is sometimes 
called CZ proper (e.g. [Ferriz-Mas, 1996 ). These differ- 
ences in the overshoot models may have some importance 
as it has been shown that flux tubes with equipartition 
field strength (~ 10** G) can be stored in all the overshoot 
layer while the field strength of 10^ gauss required in or- 
der t o explain the observe d active regions at the surface 
(e.g. Schiissler et al., 1994), can only be store d where the 
convective flux is negative (Ferriz-Mas, 199(:)- 



3. THEORY: SOME MODELS OF THE 
TACHOCLINE AND THEIR PREDICTIONS 



3.1. Purely hy drody namic models 



The first theory of th e tachocline was developed by 
Spiegel & Zahn (1992| ). They mainly address the prob- 
lem of the thickness of the tachocline, which we summa- 
rized as follows. The differential rotation ^i{9) at the 
top of the radiative interior induces a latitudinal temper- 
ature gradient and therefore a meridional circulation that 
would, without any stress acting on it, spread towards the 
interior in a thermal time-scale leading to a non rigid rota- 
tion rate in the interior and a very thick tachocline for the 
present Sun. We therefore need to invoke some stresses 
acting in the tachocline in order to prevent its progression 
towards the interior Using a purely (i.e. non-magnetized) 
hydrodynamic m odel where the tachoclin e is treated as 
a boundary layer, Spiegel & Zahn (1992 ) suggest that a 
strong enough anisotropic turbulence could lead to a thin 
tachocline. They thus obtained a relation between the 
thickness of the tachocline w and the horizontal turbu- 
lent viscosity coefficient assumed to be much higher 
than its vertical counterpart. 



8.34 10^ 



\ w / 



(1) 



Nevertheless there remain two major questions with this 
approach. First, the latitudinal rotation profile inside the 
tachocline is likely to be linearly stable to 2D distur- 
bances acco rding to Rayleigh's criterion (C harbonneau 
et al., 1999b) and so it is probably not the process lead- 
ing to the anisotropic turbulence assumed in the model. 
A solution to this problem ma y be found in a recent work 
of pikpati & Oilman (2000a ) showing that, if we allow 
radial deformations of the layer and use the so-called 
shallow water model as a first approximation of the full 
3D problem, then the tachocline is found unstable in the 
slightly sub-adiabatic overshoot zone even in the purely 
hydrodynamic case (Fig. ||). The second difficulty with 
this hydrodynamic model is that, even if it leads to a lati- 
tudinal independent rotation at the base of the tachocline, 
it cannot explain the fact that there is also no significant 
radial gradient in the interior. 
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Figure 2. Stability diagram from the shallow 
model ( Dikpati & Gilman, 2000J^ ). The variation 
G ~ 10'^ I V — V ad\ with radius is given by 
line for schematic solar model. The stability 
the tachocline thickness and location have been 



from observations { Charbonneau et ai, 1999c 



overshoot layer has arbitrarily been assignee 
ness of O.O15i?0 and it is slightly sub-. 
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3.2. Magnetized models 



Because it is very likely that the tachocline if magne- 
tized, another mechanism that can lead to anisot ropic tur- 
bulence in the tachocli ne is magnetic instab ili y. It has 
effectively been sho wn (|Gilman & Fox, 1997 



)99|_Dik- 



pati & Gilman, 1999) that there exists a joint iistability 
between latitudinal differential rotation and toro dal mag- 
netic fields. But two different kind of magnetic fields are 
usually invoked in the solar interior, namely p imordial 
fields and dynamo generated fields. As we sh; 11 see, in 
both cases theories relate their properties to the shape of 
the tachocline. 

(i) Dynamo generated magnetic fields. In viitually all 
dynamo theories (overshoot layer, int erface or Babcok- 



Leighton dynamos, see the review of [Petrova} 



the shape and especially the thickness of the 
and the overshoot layer are key issues since 
mine the strength of the magnetic field that can 
and the process by which it is transformed. 
MHD version of their shallow water model 



Gilman, 2000b) demonstrate that the presence 
netic field in the tachocline would induce a prol; 
the tachocline otherwise oblate. This can also 
by helioseismology by measuring the central 
the radial shear at different latitudes. Such 
ment of the prolateness would allow us to est 
strength of the magnetic field depending on itfe geome- 
try and localization inside the tachocline. The shape of 
the tachocline as obtained from this model is sHown as a 
function of the toroidal field strength on Fig 3 

(ii) Primordial magnetic fields. Independent y of any 
magnetic field that may be generated in the CZ, 3ough & 
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argue that we must also have a magnetic 
adiative interior in order to explain the uni- 
lo radial shear) rotation rate observed there, 
ale poloidal magnetic field would also con- 
in a thin layer as a result of a balance be- 
diffusion of the magnetic field and down- 
by the thermally driven tachocline circu- 
' obtain the following relation between the 
the tachocline and the strength of the inter- 
t field: 



Such intern d field is also required in the magnetic mod- 
;er & Kitchatinov (1997|) and M acGregor & 
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CharbonnccU (1999) where the tachocline is identified 
with an MUD boundary layer located in the radiative 
interior Assuming no magnetic coupling at the core- 
en yelope ini erface and tha t advective effects (as modeled 
by [Spiegel ^ Zahn (1992| )) dominate the viscous effects 



to prevent tl le inward spreading of the layer, M acGregor 
& Charbonileau (1999) obtain: 



3 IQ-** f — ' ^ 



(Gauss). 



(3) 



This work also suggests that the primordial magnetic field 
ly to be decoupled from the tachocline i.e. 
coHtained within the radiative interior. The in- 
bqtween a primordial time-independent mag- 
would extend into the CZ and the time- 
dynamo-generated magnetic field had never- 
investigated by ^oruta (1996 ) suggesting 
mtleraction would lead to an inversion and an 
of the primordial field by an amount that 
in on the thickness of the tachocline. On the 
if the dynamo magnetic field is generated in 
been shown to have only a weak influence 
dynatnics of both the tachocline and the radiative 



h IS 1 



(Huieirat& Yorke, 1998; jaraud, 1999) 
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4. HELIOSEISMOLOGY: TOOLS AND METHODS 
FOR INFERRING THE TACHOCLINE 
PARAMETERS 



One of the observational goals of helioseismology is to 
localize precisely the tachocline with respect to the others 
described in Sect. ^. In that perspective, we need to define 
clearly what we understand by tachocline parameters. 



4. 1 . The tachocline parameters: definition 



Following [Cosovichev (1996) we parameterize the rota 



tion rate between 0.4 and 0.8 solar radius by an error 
function of the form: 



n(r,0)^ilo^ |1 



-erf 



0.5 w{e) 



(4) 



where ilo represents the rigid rotation rate in the radia- 
tive interior, fli{6) the rotation rate at the top of the 
tachocline, and rc{0),w{0) the central position and the 
width of the tachocline at the colatitude 9. In order to 
better take into account a potential trend in the CZ an- 
other parameter(Q;) is sometimes fitted b y adding a term 
a{r - 0.7) to Eg. |4| (|Antia et al., 1998[: C orbard et al., 
1999). An important point to notice is that, with this pa- 
rameterization, the width w{9) corresponds to a change in 
the rotation rate of 85% of the jur np 17 1 (6*) — iln. O ther 



parameterizations have been used ( A.ntia et al., 1998 ) but 
we can easily compare results using this later definition. 



4.2. The inverse problem: principle and difficulties 

Solar oscillation modes are identified by three integers: 
the spherical harmonics degree I and azimuthal order to, 
and the radial order n. The frequency splitting Svnim 
induced by the rotation il(r, /x) can be calculated by: 



(5) 



where, r is the fractional solar radius, /i = cos{9) and 
Knim are kernels calculated from a standard solar model. 
Inferring the internal rotation rate from observed split- 
tings therefore requires inverting this integral equation. 
In order to achieve this, two classes of methods, both 
linear, are usually used. The first approach i s a Least- 
Squares (LS) method (e.g. [Corbard et al., 1997 ) by which 
we try to find the rotation profile that minimizes the sum 
of the square of the differences between observed and 
predicted splittings, weighted by the observational errors. 
The second app roach, called Optimally Loc alized Aver- 
ages (OLA, e.g. Pijpers & Thompson, 1992), consists in 
trying to reach locally, at (rg, I-Iq), the best resolution by 
building the appropriate linear combination of observa- 
tions < n{ro,i'o) >^J2cnir,i{ro,m)Snirn- From Eq. | 
this is equivalent to taking a linear combination of the 
kernels A'„;,„. The resulting kernel 



is called averaging kernel because it follows from the 
previous equations that the quantity < il(ro,/io) > is 
an average of the rotation rate in a domain defined by 

K{ro,iJ,o,r,fi): 



<fl{ro,lia) >= / K{rQ,fj.o,r,fj.)^{r,fi)drdfi 
Jo Jo 

(7) 

The difficulties arise from the fact that Eq. |] is an ill- 
posed problem with no unique solution. In the global 
(LS) methods, we therefore need to introduce some a- 
priori knowledge on the rotation profile in order to regu- 
larize the solution and avoid strong oscillations. Never- 
theless, this Regularized Least Squares (RLS) approach 
prevents us from recovering accurately sharp gradients as 
those expected in the tachocline because they do not con- 
form to the global smoothness a-priori introduced. 

For local (OLA) methods, the limitation in resolution 
comes from the error propagation. A quantitative mea- 
sure of the resolution in the radial direction is obtained 
by fitting the averaging kernel by a Gaussian profile cen- 
tered at tq with a standard deviation ctj.- We then define 
the radial resolution by ~ y^So^. The reason for the 
choice of this definition will become more clear in the 
next section but it turns out that from actual observations 
the tachocline is not resolved by linear inversions when 
limiting error propagation to a reasonable amount. 

Therefore the information about the thickness of a sharp 
gradients of the rotation profile is not directly readable 
from the solutions obtained by classic inversions. How- 
ever three different approaches have been developed in 
order to overcome this difficulties. 

(i) Forward analysis. The general shape of the rotation 
rate is assumed to be known and the tachocline is parame- 
terized with few parameters (as the erf function of Eq. ^ 
whic h are adjusted to fit the data by calibration (B asu, 
1997) or by min imization using m ethods such as simu- 
l ated annealing ( ^ntia et al., 1998 ) or genetic algorithm 
(phai-bonneauet al., 1999a|). 



(ii) Deconvolution. This method uses our knowledge of 
the resolution kernel in linear methods in order to reach 
a better estimate of the tachocline width. The basic idea 
is to approximate Eq. |7|by a convolution equation (C har- 
bonneau et al., 1999a; [Corbard et al., 1998| ). Then, if the 
tachocline profile after inversion is approximated by an 
erf function of width w, the 'true width' Wc can be ob- 
tained by: 

(8) 



A2 



where A,, is the radial resolution at the center of the 



tachocline as defined in Sec. 4.2 



(iii) Non linear or ada ptive regularization. T his approach 
has been su ggested by [Corbard et al. (1998 ) and fully de- 
veloped by Corbard et al. (1999). In brief, we construct 



\cnhn{rQ,^lo)Knhn{r,tj) (6) 



an iterative process which will adjust locally the smooth- 
ness term as a function of the gradient amplitude found 
at the previous step. This allow us to keep the well con- 
strained sharp gradient zones while still regularizing else- 
where. 
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Table 1. Tachocline width and position as inferred from different authors, data and methods. The upper table gives values 
over an average in latitude (the weighting function being (1 — /i^)(5/i^ ~ 1)) while the lower table is for fits at the equator 
GONG month 4 starts 8/23/95. The end of GONG months 7,10 and 14 are respectively 1/13/96, 4/30/96, 9/21/96. LOWL 
data have been collected between 2/26/94 and 2/25/96. 



Authors 



TJR 



■o 



■o 



Data 



Method 



Kosovichev(1996) 



Basu (1997) 



[Charbonneau et al. (1999a|) 



0.692 ±0.005 
0.7066 ±0.0047 
0.7034 ± 0.0056 
0.7048 ± 0.0039 
0.705 ±0.002 



0.09 ±0.04 
0.0412 ±0.0260 
0.0490 ±0.0245 
0.0514 ±0.0177 
0.053 ±0.015 



BBSO 86,88-90 
BBSO 86,88-90 
GONG 4-7 
GONG 4-10 
LOWL 94-96 



FM 

FM (calibration) 
FM (calibration) 
FM (calibration) 
Genetic FM 



Weighted average 



0.704 ±0.002 



0.052 ±0.010 



[Charbonneau et al. (1999a[ ) 



Corbardet al. (1998 ) 



L^orbard et al. (m% 



Antia et al. (1998) 



0.689 ±0.006 
0.691 ±0.007 
0.695 ±0.005 
0.691 ±0.004 
0.6851 ±0.0077 
0.6843 ±0.0112 



0.01 ±0.03 
0.07 ±0.03 
0.05 ±0.03 
0.01 ±0.03 
0.0230 ±0.0407 
0.0098 ±0.0093 



LOWL 94-96 
LOWL 94-96 
LOWL 94-96 
LOWL 94-96 
GONG 4-14 
GONG 4-14 



OLA + Deconvolution 
RLS + Deconvolution 
RLS + Deconvolution 

Adaptive Regularization 
FM (calibration) 

FM (simulated annealing) 



Weighted average 



0.691 ±0.003 



0.034 ±0.014 



RESULTS FOR THE EQUATORIAL PLANE AND 
FOR SPHERICAL AVERAGE 



All modes have amplitude near the equator and there- 
fore we expect to have more chance to be able to infer 
tachocline parameters accurately there. The first results 
on the tachocline parameters were obtained (Tab. |l]) ei- 
ther at the equator or on an average over all latitudes. 
The equatorial results leads to an interval [0.674,0.708] 
or, taking into account uncertainties on both position and 
with, [0.665, 0.717] for the location of the tachocline. In 
the case of the latitudinally averaged results, we obtain 
respectively [0.677, 0.730] and [0.671, 0.736]. These first 
results show that the equatorial tachocline is centered sig- 
nificantly is^ O.O2i?0) below the base of the CZ and 
that it is very likely that all the strong radial gradient 
of angular velocity is located below in the equatorial 
plane. On the other hand, the latitudinally averaged re- 
sults suggest that these characteristics are not maintained 
at all latitudes leading, on average, to a thicker tachocline 
centered only around O.Oli?© below Tcz- Therefore, it 
seems that, in a latitude-average, the upper third of the 
tachocline is located in the nearly adiabatically strati- 
fied zone above Tcz- Nevertheless this results must be 
taken with caution because the decreasing radial reso- 
lution with latitude tends naturally to show a thicker 
tachocline at high latitudes thereby influencing the aver- 
age thickness and position. The parameters obtained this 
way should therefore be regarded as upper limits for the 
tachocline properties. 

The latitude-average width is however probably the best 
estimate to use in the formulae obtained from the various 
theories of Sect. ^ which assume also a spherical sym- 
metry. In the hydrodynamic theory, Eq. [l| with a with of 
O.O5i?0 leads to an horizontal turbulent viscosity coeffi- 
cient of about 3 lO^cm^ s^^, several order of magnitude 
higher than the microscopic value. In the MHD theories 
the same value of the width leads to a primordial mag- 
netic field strength of [Sol ^ lO^^'Gauss for both Eq. g 
and H This suggests that even a weak magnetic field is 



enough to keep the radiative interior rotating rigidly and 
to confine the radial shear to a thin layer compatible with 
observations. By analyzing the even splitting coefficients 
from GONG observations (see Sect, ^j] below), B asu 
(1997) set an upper limit of 0.3MG for a field located 
at the base of the CZ. From Eq. ^ this would correspond 
to a width of O.OO45i?0 which cannot be excluded from 
observations . By shearing a poloidal field of 10 ^ ''Gauss 
the theory of MacGregor & Charbonneau (1999 ) predicts 
the generation of a toroidal field of ^ O.IMG which is 
also compatible with the 0.3MG limit. 



6. VARIATIONS WITH LATITUDE? WITH TIME? 



Theories predict variation with latitude of the shape of 
the tachocline. We have seen for example that, in the 
shallow-water model, the presence of a magnetic field 
would induce a prolat e tachocline otherwise oblate. Both 
the magnetic model of tludiger & Kitchatinov (1997) and 
the 2D hydrodynamic model of Spiegel & Zahn (1992 ) 
predict a thicker tacho cline at the pole than at the equator 
(see also Elliott, 19971). In the model of G ough & Mcln- 
tyre (1998) the tachocline is thicker where the horizontal 
component of the magnetic field in the radiative interior 
is small. The poles should be therefore thicker but that 
may also be true for other latitudes. This means that the 
shape of the tachocline could also be a diagnostic for the 
geometry of the magnetic field. Some attempts have been 
made to infer latitudina l variations in the t achocline pro- 
file from o bservations (|Antia et al., 1998|; C harbonneau 
et al., 1999a), but no compelling evidence has been found 
yet. Even if, as shown in the previous section, the re- 
sults tend to argue in favor of a prolate tachocline thicker 
at high latitudes than at the pole, this may very well be 
due to the limits in resolution associated with the datasets 
used. 



Time variations are also expected mainly because of the 
changes of the magnetic field strength and geometry dur- 
ing the solar cycle. Some evolution of the large scale 
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100 1000 10000 

v/I (/iHz) 

Figure 4. Distance in frequency between neighboring i's 
at constant n as a function ofv/L. The dashed lines show 
the frequencies corresponding to the first three sidelobes. 
The frequencies used here are an average over six inde- 
pendent data sets. 



circulation inside the layer or the action of internal grav- 
ity waves may also induce temporal changes or oscilla- 
tions in the tachocline parameters. Small amplitude os- 
cillations, with a period of 1.3 year, have recently been 
found in the rotation profile clo se to the boundaries of the 
tachocline ( ^iowe et al., 2000 ). These results are based 
on both GONG and MDI data but have not been c on- 
firmed by independent analysis (A.ntia & Basu, 2000). If 



these oscillatio ns are real, the 1.3 year period remains un- 
explained (see Gough, 2000). Because of the observed 
magnetic cycle, we rather expect the tachocline variations 
to be on the time scale of the solar cycle (11 years). A 
first attempt at detecting long term temporal variations in 
the tachocline have been made by Basu & Schou (2000[ ) 
using 1 1 sets of MDI observations covering 72 days each 
between July 1996 and April 1998. They found no clear 
changes in the width or position of the tachocline but gave 
some hints that the position may move slightly outwards 
with increased activity. In the following we use the obser- 
vations collected by LOWL instrument during the 6 years 
of the ascending phase of the solar cycle (1994-2000) in 
order to investigate the possible cycle related changes in 
the tachocline. But lets first introduce this new dataset. 



6.1. Six years of LOWL observations: a new dataset 



The LOWL instrument ( lomczyk et al. , 1 995 ), based on a 
Potassium magneto-optical filter, has been collecting re- 
solved Doppler observations for more than six years. The 
spectra were processed usin g the LOWL pipeline, which 
has been recently improved (Jimenez-Reyes, 2001 ). An- 
nual time series for degrees from £ = up to 99 have 
been created by using spherical harmonic masks and a 
Fast Fourier Transform has been applied to each one of 
the time series. The average duty cycle over one year of 
observations was around 20%. 

A detailed description of the fit ting method can be found 
in Jimenez-Reyes et al. (2001). In brief, we have used 



the general expression of the likelihood function, assum- 
ing that the statistics of the real and imaginary parts of 
the Fourier transform follows a mu lti-normal dis tribution, 
described by a co variance matrix (ISchou, 1992|; Appour- 
chaux et al., 1998). This model is required due to the fact 




[(3>'/Sl).-11.57]/r 



Figure 5. Average of ai over six years. The solid line de- 
notes the best fit of a Gaussian profile plus a background 
given by a straight line. 



that the spherical harmonics are not orthogonal over the 
observed hemisphere. This introduces a natural leakage 
between modes with different £ and m. 

The observed splitting between m-components is given 
in terms of a-coefficients by: 



E 



Vntm = Vnt+ } ^ aj{n,£)P, (m) 



(9) 



where P^{m) are orthogonal polynomials normalized 



such that Pf{£)=e ( |Schou et al., 1994 App. A). The 



sum over the odd a-coefficients gives the rotational split 
ting Snim (Eq. |^, while the even terms are mainly due to 
magnetic fields and second order effects of the rotation. 
An analysis of the central frequency Vni and the even 
a-coefficients has been carr ied out using these data by 
( Fimenez-Reyes et al., 2001). These parameters present 
significant variations very well correlated with the solar 
cycle. 

One of the main source of systematic errors in the fit- 
ting procedure is the leakage, which affects mainly the 
a-coefficients estimates. It is particular important when 
the distance in frequency between neighboring £'s at con- 
stant n gets smaller Following the asymptotic expres- 
sion given by e.g. Deubner & Gough (1984) for the 



frequencies at intermediate and high degrees, we obtain: 
~ Therefore, the distance between neighbor- 
ing ^'s is proportional to i^/£ which, in turn, is related to 
the inner turning point radius rt through the local sound 
speed c{rt) by: 27rrt = c{rt){L / v). This property is 
illustrated on Fig. ^ Moreover, in the case of observa- 
tions collected fromjust one site (as LOWL), sidelobes at 
1 1.57/iHz from the main peak will appear in the Fourier 
Transform due to the modulation of one day introduced 
in the signal. The leakage between modes will therefore 
occur at lower turning point. Using Fig. ^ we can pre- 
dict where to expect complications in the fitting. The 
dashed lines denote the position in frequency of the first 
three sidelobes. When v/L k, 31/iHz the first sidelobe 
from l±l will cross the position in frequency of the tar- 
get mode with degree I. The amplitude of the second 
and third sidelobes are very small will not be considered. 
From the equation above and a standard sound speed pro- 
file, we can deduce that v/L 31/iHz corresponds to 
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Figure 7. Difference between the central position infered 
at 60° and at the equator, as a function of time. The upper 
plot corresponds to the results of Fig. ^ while the lower 
plot corresponds to inversion using stronger regulariza- 
tion. 



Figure 6. Inferred tachocline parameters as a function 
of time from LOWL observations. The full and dashed 
lines correspond respectively to fit at the equator and at 
60°. For the width (upper left), the lower dashed line is 
obtained after deconvolution using Eq. The scale on 
the right of the lower right panel corresponds to the fit at 
60° (dashed line). 



Tt = O.9i?0 which is effectively the depth where the 
previous analysis of LOWL data s howed an artificial in- 
crease of the solar rotation (e.g. Corbard et al., 1997 ; 
Chaplin et al., 1999 ). Nevertheless, two different peaks 
very close in frequency can be overlapping or not depend- 
ing on their linewidth F. Therefore, the a-coefficients 
should be plotted against [((5i^/W)„ - 11.57]r-i (Fig.|). 
An important feature or hump appear centered approx- 
imately at zero. Its amplitude is maximum for ai 
2nHz). The same analysis was carried out year by year, 
showing that this feature is found similar for all time se- 
ries. Thus, we decided to remove this systematic error by 
fitting the bump with a Gaussian profile in the average 
over six years and by subtracting afterwards this profile 
from every yearly fit. The even a-coefficients do not show 
similar humps and have therefore not been changed. 



6.2. Results and Discussion 



The splitting coefficients obtained for each year have 
been invert ed using a 2D RLS inversion code (C orbard 
et al., 1997). The rotation profile has then been fitted at 
the equator and 60° by an erf function as given by Eq. ^. 
The averaging kernels at the center of the tachocline have 
been fitted by a gaussian in order to estimate the radial 
resolution. The results are shown in Fig. ^. The radial 
resolution obtained is about O.12i?0 at 60° and O.lOi?© 
at the equator The angular resolution achieved is related 
to the number of a-coefficients in Eq. ^ Our analysis 
includes up to 9 coefficients leading to a maximum an- 
gular resolution of about 20° at the equator. Because the 
radial shear doesn't exist at about 30° and because the 
angular resolution decrease at high latitudes, we limit our 
analysis to two latitudes: the equator and 60°. There- 



fore, in the following, prolateness refers to the difference 
between the central position of the layer at these two lat- 
itudes. 

Because the inferred width at the equator is always lower 
than the resolution we cannot use our simple model 
for deconvolving. This indicates that the width of the 
tachocline at the equator is probably lower than the local 
spacing of the grid used for the inversion i.e. O.O2i?0. 
The same happens at 60° for 1995, 1998 and 1999. 
There is therefore no strong evidence of a systematically 
thicker tachocline at 60°. Moreover the errors reported 
on the plot are formal errors as obtained from the fits 
but Monte Carlo simulations have been carried out which 
suggest that the uncertainties on the inferred widths af- 
ter the whole inversion process are between ±0.02 and 
±O.O3i?0. 

The center of the tachocline is always found deeper at 
the equator than at 60°. The variation of this prolateness 
with time is shown in Fig. ^ The maximum of prolate- 
ness found is about O.O2i?0 and no prolateness is found 
the first and last years. These differences and their fluctu- 
ations are nevertheless very small and are also very sen- 
sitive to the inversion parameters chosen and especially 
the amount of regularization used. The lower panel of 
Fig. illustrates this point by showing that, with a more 
regularized inversion, the maximum prolateness is about 
O.O3i?0 and a minimum is no longer found for 1994. 



Generally speaking, we do not find from this analysis ev- 
idence of any general trend or significant oscillation in 
the tachocline parameters during the ascending phase of 
the actual solar cycle. In particular we do not find an out- 
ward trend for the central position as suggested by B asu 
& Schou (2000). Nevertheless, a prolateness of the layer 
is observed every year which allow us to, at least, exclude 
an oblate tachocline. The amount of prolateness found 
between the equator and 60° is however of the same or- 
der of magnitude than the uncertainties on the width of 
the layer. We can therefore only set an upper limit for 
the prolateness which is around 0.03i ?f7i. This is in good 
agreement with previous estir nates of Antia et al. (1998 ) 
(0.004<A^c/i?o<0.023) and pharbonneau et al. (19'99a) 
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(Arc — O.O23i?0). Using the shallow-water model and 
shape curves as shown in Fig. 3, this would correspond 
to a toroidal magnetic field strength of about O.IMG if 
it is located in the overshoot layer or about 0.6MG if it 
is located in the radiative interior. If the toroidal field is 
concentrated in bands migrating towards the equator dur- 
ing the ascending phase of the cycle, one would expect, 
from the shallow water model, a decreasing prolateness 
for the period of LOWL observations. This general trend 
cannot be excluded but is not observed from our prelimi- 
nary analysis of LOWL data. 
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